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A broken symmetry ground state without any magnetic moments has been calculated by means 
of local-density-approximation to density functional theory plus a local exchange term, the so-called 
LDA+[/ approach, for URu2Si2. The solution is analysed in terms of a multipole tensor expansion 
of the itinerant density matrix and is found to be a non-trivial magnetic multipole. Analysis and 
further calculations show that this type of multipole enters naturally in time reversal breaking 
in presence of large effective spin-orbit coupling and co-exists with magnetic moments for most 
magnetic actinides. 



The magnetism involving the 5/ states of the actinides 
shows many exotic behaviours. One of most enigmatic 
is the tiny magnetic moment [TJ [5] observed in tetrag- 
onal URu2Si2 below a transition temperature of 17 K. 
It is established that the true ordering parameter driv- 
ing the second order transition is hidden, i.e. not ob- 
servable by standard techniques. A huge experimen- 
tal effort has been made to uncover the true nature 
of this hidden-order. This has led to a large progress 
in determining the relevant phase diagram of URu2Si2, 
but the order stay concealed. Under pressure the or- 
dering temperature increases slightly, but more signifi- 
cantly for pressures above 0.5 GPa there is a phase tran- 
sition to an anti-ferromagnetic (AF) state with magnetic 
moments of 0.3 /j.b[3]- The critical pressure increases 
with temperature and seems to reach a bicritical point 
around 1.3 GPa. The two phases HO and AF have some 
properties in common. They both show large anoma- 
lies in thermodynamical0] and transport quantities |S] at 
the critical temperature, signalling a large Fermi surface 
nesting. In fact these quantities^ as well as the Fermi 
surface geometry determined by de Haas-van Alphen [B] 
show only smooth variations across the phase transition. 
Since there is a sharp transition between the two phases, 
and since the tiny moments observed in the HO phase 
are argued to be caused by sample dependent extrinsic 
effects, the order parameter in the HO can in principle be 
of completely different symmetry than in the AF phase. 
Especially, the HO might be time reversal (TR) symmet- 
ric although the AF clearly is not. On the other hand, 
the inflection point observed in the induced moment of an 
external magnetic field, indicates that there is an "adia- 
batic continuity" between the two phases and hence that 
HO is TR oddilj. In summary, several crucial aspects of 
the HO are still under debate: whether the HO is of even 
or odd TR symmetry, whether it is of localised or itin- 
erant character, and whether the HO and AF phases are 
related or of different origin. Many explanations, more 
or less exotic, have been suggested for this "hidden or- 
der" (HO), some of localised nature e.g. quadrupoles [8j 
or octupoles^, and some of itinerant nature as e.g. spin 
nematics[TU], orbital currents [TT]. helicity order |T2|, or 



fluctuating moments [13], but none of them have been 
consistent with all experimental observation in a satis- 
factory way. 

In this Letter we find that a non-trivial magnetic tri- 
akontadipole (rank five) moments constitue the HO pa- 
rameter in URu2Si2. In fact we observe that these mul- 
tipole moments play an important role in most mag- 
netically ordered materials with a large effective spin- 
orbit coupling (SO). Our calculations which treat the 5/ 
states as itinerant with large Coulomb interaction show 
that what makes URu2Si2 unique is the very large tri- 
akontadipole moment which forces both spin and orbital 
dipole moments to vanish. In this material the Fermi 
surface nesting is found to be important for the stabi- 
lization of the calculated staggered multipole moments. 
With decreased in-plane lattice constant the dipolcs in- 
crease giving rise to ordinary magnetic moments in ac- 
cordance with pressure experiments [3] as the magnitude 
of the magnetic multipole decreases slightly. 

The driving mechanism for magnetic ordering in con- 
ventional magnets based on id transitional metal ele- 
ments, is known to be the Stoner-like exchange, which 
is for instance included in local density approximation 
(LDA) to the spin density functional theory (SDFT) 
approach [TJ. This gives in most cases a very good 
description of magnetic moments and magnetic order- 
ing, whether ferromagnetic, anti-ferromagnetic or non- 
commensurate spin density waves. When relativistic ef- 
fects, especially the SO coupling, start to play an im- 
portant role as in the actinides, the Stoner-like spin- 
polarisation is not sufficient anymore. It tends to over- 
estimate the spin contribution to the magnetic moments 
while simultaneously it drastically under-estimates the 
orbital contribution. This is known to be remedied by 
adding a more general exchange terms, that includes for 
instance an orbital polarisation term|15[ I16j . A gen- 
eral form of exchange interactions within an atomic open 
shell is the screened Hartree-Fock term. This term is 
fused with the conventional SDFT formalism in the so- 
called LDA+C/ method [T?! [18]. Recently we have shown 
that with this type of exchange, the SO coupling can 
be strongly enhanced, i.e. the dominating tensor mo- 
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ment contribution comes from the multipole tensor w^^*^. 
In the case of i5-Pu this was shown to quench the spin- 
polarisation leading to a TR even state [19] in accordance 
with experiments. The analysis of the ground state was 
achieved by a general multi-pole expansion of the ex- 
change energy, Ex ~ J2kp7- ^kpr'^^^^ ■ w'^p'', where w'^p'' 
is the multipole tensor moment for the channel kpr and 
Kkpr is the corresponding exchange parameter These 
tensors correspond to density multipoles for even 
k even and current multipoles for odd. The index p = Q 
describes the charge density and current while p = 1 de- 
scribes the spin density and curent. The orbital {k) and 
spin {p) indices are coupled to the index r. Hence w^^'^ 
corresponds to a scalar product of spin and orbital cur- 
rent s ■ £, while e.g. w°^^ and w^'^^ are proportional to 
spin and orbital moments. 

In this study we apply the LDA-I-J7 approach and its 
multipole decomposition to a few different magnetic ura- 
nium systems, three supposedly normal systems; UAs, 
US and USb in the NaCl-structure, and URu2Si2 with 
its enigmatic HO in the tetragonal ThCr2Si2-structure. 
The calculations are set up according to the observed 
crystal structure [Ij [20] as well as the observed magnetic 
ordering wave vectors; a single q — (00^), q — (000) 
and a triple q = (00|), respectively for the three "nor- 
mal" compounds. The calculations where done in the 
APW-|-to[lT] method together with the LDA+U scheme 
as implemented [m [22] in the code Elk [23]. The basis 
cut off Rl]Gms^x, where i?u is the muffin tin radius of ura- 
nium, was set to 9.5. The integration over the BZ zone 
was performed with a 12 x 12 x 12 /c-points for the fee cell 
of the ferromagnets UAs and US, with 8x8x8 /c-points 
for the simple cubic cell used in the triple q ordering of 
USb , and with 18 x 18 x 10 for the simple tetragonal cell 
of the single q ordering of URu2Si2. 

The calculated magnetic moments are displayed in 
Fig. [l] as a function of the parameter U. We notice 
that we have a good agreement with the experimental 
moments [^ for U in the range of 0.6-1.0 eV, which is 
the range of reasonable values for the 5/ states of U. 
These moments are also of the same magnitude as in 
other "beyond-LDA" calculations [13 IHj. The calculated 
spin and orbital moments (SM and OM) show slightly 
different behaviour for the different uranium compounds, 
for UAs the OM increases monotonously with U and the 
SM is constant while for US and USb the OM has a max- 
imum value around 0.6 eV and the SM monotonously de- 
crease. The last behaviour is somewhat counter-intuitive 
- the SM decrease with increased exchange energy. Here 
we should note that we have adopted a strategy where 
the different Slater parameters are screened equivalently 
[55[[2S]. This has advantage that the limit of vanishing 
U corresponds to an ordinary LDA calculation. 

In order to understand why the spin-polarisation de- 
creases with U, we display the leading tensor contribu- 
tions to the exchange energy [THj [55] as a function of the 
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FIG. 1: At top: The calculated spin (black), orbital (red) and 
total (green) magnetic 5/ moments as a function of U, for var- 
ious metallic and magnetic uranium compounds. U = corre- 
spond to the LDA-limit. The experimental magnetic moment 
is given by the dashed line. At bottom: The contributions to 
the exchange energy from the relevant polarisation channels 
kpr as a function of U. 

screened [/-parameter in Fig. [T] We can see that the 
spin-orbit coupling like multipole w^^" has a large con- 
tribution to the total exchange energy. This leads in 
turn to that the spin-polarisation decrease in importance 
compared to the LDA-limit, U = 0. Although the finite 
U give rise to an orbital polarisation, the most striking 
effect is that the a non-trivial high order multi-pole dom- 
inates, the w^-'^^ tensor, and that it plays a large role for 
all the uranium compounds. As given by its kpr indices 
this contribution arises from a multipole of the spin mag- 
netization density. 

This observation of large higher order multi-poles, nat- 
urally leads to the question: "how relevant are they 
for URu2Si2 with its tiny magnetic moments and hid- 
den order parameter?" In Fig. [T] we have also display 
the calculated moments and decomposed exchange ener- 
gies for URu2Si2, set up with an antiferromagnetic order 
q = (001). Here we see that while the overall tendency is 
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FIG. 2: In the upper part, the spin, orbital and total mag- 
netic moments, are shown as a function of variation in the 
in-plane lattice constant a with respect to the experimental 
one ao, calculated with U = 0.9 eV. In the lower panel the 
non-vanishing contributions to the exchange energy are dis- 
played. 



in accordance with what we found for the other actinide 
compounds above, there are two details worth pointing 
out. Firstly, the polarisation of the 615-channels are even 
more prominent and, secondly, at large U (1.2 eV) the 
SM switches sign and becomes parallel to the OM and 
both moments become very small. These properties sig- 
nals some anomaly, and we proceed to study this phase 
under pressure. From careful analysis of the results of ex- 
periments under uni-axial stress, one has concluded that 
the main variation in the pressure experiments arise from 
the contraction of the tetragonal a-axis[26 . Therefore 
we have performed calculations with varying lattice con- 
stant a, to mimic the effect of pressure. The results for 
U — 0.9 eV are shown in Fig. [2] where we see a dramatic 
effect on the magnetic moments from small variations 
in lattice constant a away from a critical value, slightly 
(1.4 %) larger than the experimental value aq. At this 
critical value both the SM and OM vanish, while simul- 
taneously the 615 contribution to the exchange energy 
almost diverges. In fact all the vector contributions; Oil, 
101 and 211, go to zero at this point. The magnetic mo- 
ments only disappear for a narrow range of the in-plane 
lattice constant in the present study where only this pa- 
rameter is allowed to vary. Future studies have to clarify 
whether this range increases when proper total energy 
optimisations, including proper Si site relaxation as well 
as variation of the c-axis, are taken into account. 

These results are very sensitive to the Brillouin zone 



integrations, which indicates that the Fermi surface 
nesting |13j of the un-polarised state play a major role in 
this instability. This nesting will lead to a large degener- 
acy that can be split by TR breaking polarisations. Since 
in the present case the effective spin-orbit coupling is very 
large this TR breaking cannot be of usual spin-splitting. 
Instead it is more convenient to study this TR breaking 
in the j = 5/2 subset of the / shell. Then it is possible to 
see that it is order parameters of the tensor type w'^^'^'^^^) 
that occurs after the symmetry breaking, with larger con- 
tribution the larger k. This is why the 615 tensors mo- 
ments enter. Further analysis gives that there are two 
components of the 615 tensor that are allowed by symme- 
try. Both are large but it is the w'^i = {wf-^ + w'l^i) / V2 
component that drives the rapid variation in exchange 
energy as function of the in-plane lattice constant, while 
Wq^^ are almost constant. The spatial variation of the 
magnetization density arising from these triakontadipole 
components are displayed in Fig. [3j From this plot we 
can note that it is closely related to the intra-atomic non- 
coUinear spin density that always arise in the presence of 
spin-orbit coupling [27] . Since the exchange coupling en- 
hances the effective spin-orbit coupling substantially, the 
resulting 615-polarisation becomes more important than 
the ordinary spin-polarisation Oil. For all cases except 
UAs its exchange energy contribution bypass the one of 
spin-polarisation for values of U around 1 eV. This casts, 
together with studies on other magnetic actinides|22j. 
some new light on the magnetism of uranium systems 
and actinides in general. 

From Fig. |2] we see that for smaller lattice constant, 
a, antiferro-magnetic moments rapidly occur. For low 
pressures the SM takes a value of about 0.5 /ie, while 
the OM is larger and of opposite sign. This state is in a 
sense anomalous too since the SM is now anti-parallel to 
the Wq^^ moment, while they are parallel in the limit of 
weak spin-orbit coupling. However as is clear from our 
calculations, this AF order is of minor importance since 
it is still the HO, the 615 multi-pole, that dominates the 
exchange energy. Although the fact that the AF state 
stabilised under pressure is closely connected to the HO 
phase in line with the concept of adiabatic continuity [7], 
there are two different solutions in the calculation and 
the transition from one to the other sensitively depends 
on the value of U in our calculations. For values of U 
above 1.2 eV the transition is clearly of second order, 
while for values below 0.8 eV it is of first order. In be- 
tween these values the numerical accuracy is not enough 
to safely determine the order. However, in order to prop- 
erly describe this phase transition, accompanying lattice 
relaxations have again to be fully taken into account, 
which is out of scope of the present study. 

Can the 615-tensor order parameter be observed at this 
critical lattice point where dipole tensors vanish? Well, 
due to its high rank it is indeed a well hidden order pa- 
rameter, but since there is a magnetization density asso- 
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FIG. 3: The angular variation of the direction of the spin 
density for the two non-vanishing components of the triakon- 
tadipole, Wq^^ and w;±4*. The green (red) colour indicates the 
regions where the spin axis is outward (inward) normal. 



larisation, are forced to vanish. Under pressure, the 615 
polarisation decreases and dipole polarisations are again 
allowed. In experiments this manifests as an apparent 
dipole AF order, although we find in our calculations 
that the HO candidate w^^^ still dominates the physics 
of the 5/ shell. Our results imply that there is always an 
hidden order in the magnetic actinides, it is only when 
it forces the dipoles to vanish as in URu2Si2 it becomes 
perceptible. 
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ciated with it, as in Fig. [3l it will give rise to magnetic 
scattering in e.g. neutron diffraction (ND) experiments 
although the integrated moment is zero. One problem 
though is that it belongs to the same point group repre- 
sentation as any non- vanishing dipole order [25] . so very 
careful analysis on high accuracy experiments is needed 
to distinguish this pure HO case from a tiny dipole mo- 
ment case. However, if it is this triakontadipole-state 
which is the HO, it can resolve a puzzling discrepancy 
between the tiny moments observed in ND[5H1 with 
the lack of local spin splitting detected by e.g. NMR 
experiments ^30. : It is not tiny dipole moments that give 
rise to the spin-flip observed in the ND experiments, but 
the non-coUinear 615-multipoles, which by nature have 
an extremely short range stray field and hence leads to 
no hypcrfinc-fields at the probing sites. Hence, there 
would be no need to invoke the idea of inhomogeneous 
ordering 30J, i.e. that AF order exists in part of the sam- 
ple while HO exists in the rest. 

In the calculations of the magnetic uranium com- 
pounds we generally find a coexistence of polarisations 
of the 615 channel and the Oil channel. An analysis 
of the properties of the density matrix gives that a Oil 
polarisation is always permitted, and usually favourable, 
whenever full saturation of the other channels have not 
been reached. However, when both the wl^'^ and Wq^^ 
tensor components are saturated, an additional satu- 
ration of leads to that the density matrix would 
acquire negative eigen-values by any spin polarisation 
Wq^^. Hence, these two latter types of polarisation are 
competing, in accordance with our calculated results for 
URu2Si2 in Fig. [2] In this respect, the observation of a 
HO of URu2Si2 is the effect of that the polarisation of 
is nearly optimal and hence forbids the usual spin 
polarisation. 

In conclusions, we find that the triakontadipoles w^^^ 
play a major role in all magnetic light actinides. What is 
unique in the case of URu2Si2 is that its polarisation is 
so large that the usual dipole polarisations, e.g. spin po- 
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